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We compute, for the first time, the absorptive part of the massless correlator of two quark scalar 
currents in five loops. As physical applications we consider the G(cQ) corrections to the decay rate 
of the Standard Model Higgs boson into quarks, as well as the constraints on the strange quark 
mass following from QCD sum rules. 



PACS numbers: 12.38.-t 12.38.Bx 12.20.-m 

INTRODUCTION 

Within the Standard Model (SM) the scalar Higgs bo- 
son is responsible for the mechanism of the mass gener- 
ation. Particularly interesting for the observation of the 
Higgs boson with an intermediate mass Mh < 2M\y is 
the dominant decay channel into a pair of bottom quarks, 
H — > bb. The decay of the Higgs boson into a quark- 
antiquark pair (//) proceeds through its coupling to the 
corresponding quark scalar current and reads (for a re- 
view see e.g. 



r(H -> //) 



4V2tt 1 



M 2 H ), 



(1) 



with R(s) — ImII(-s — ie)/(2irs) standing for the ab- 
sorptive part of the scalar two-point correlator: 



H(Q 2 ) = (4tt) 



dxe lqx ( 



T[J%x)J%Q)]\Q). (2) 



Here Q 1 = —q and Jy = V E'/ V I / / is the scalar current for 
quarks with flavour / and mass m/, coupled to the scalar 
Higgs boson. 

The currently known fixed order perturbative predic- 
tions for R can be shortly summarized as follows 0, El 
(we have put the number of effective flavours n/ = 5) 



R=l + 5.66667 a s + 29.1467 a? + 41.7576 a 3 . 



(3) 



with a s = a s (Ms)/7T. Note that eq. © is given in the 
high energy limit with all power-suppressed terms of or- 
der rrif/Mjj and higher neglected. In fact, the full depen- 
dence on the quark mass m/ is known up to and including 
the 0(a 2 s ) contribution We will not discuss the power 
suppressed terms in the present publication. 

For Mh > 70 GeV already the term of order 
m\ /M H oq is numerically by an order of magnitude less 



than the massless 0(a 3 s ) term displayed in eq. (|3J 3]. 
Due to the large mass Mh the couplant a s is less than 
0.04 which results in a good apparent convergence of the 
perturbation series in eq. ©. 

On the other hand, for energy scales, say, of order of a 
few GeV's, relevant for QCD sum rules the higher order 
corrections to the correlator J2J are numerically quite 
important (see, e.g. 0, IE S U)- We will discuss this 
issue in some detail later in the text for an example of 
bounds for the light quark masses. 

The motivation of the present publication is fourfold. 
First, presenting the first complete 0(cvf) result for a 
massless QCD correlator shows that the new theoreti- 
cal methods used in 0, 0, indeed do deliver genuine 
QCD results in five-loop approximation. Second, the re- 
sults are important for the QCD sum rules based on a 
(pseudo)-scalar correlator (J2J; they also provide an ac- 
curate prediction of the Higgs decay rate into hadrons. 
Third, the case of the scalar correlator should be consid- 
ered a necessary preparation step before computing the 
O(o.g) contribution to the vector correlator. The impor- 
tance of the latter calculation for the precise determina- 
tion of the value of a s from the r-lepton and Z— boson 
decay rates is well-known. Last, by comparing the exact 
af result with the estimates based on various optimiza- 
tion procedures one obtains important insights into the 
quality of different approaches, confirming some and re- 
futing others. 



CALCULATION AND RESULTS 

To compute the absorptive part of n we proceed in full 
analogy with our previous calculations described in |l2l | . 
First, u sing the criterion of irreducibility of Feynman in- 
tegrals [l3j [lj] , the set of irreducible integrals involved 
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in the problem was constructed. Second, the coefficients 
multiplying these integrals were calculated as series in 
the l/D — > expansion with the help of an auxiliary in- 
tegral representation [l5| . Third, the exact answer, i.e. 
a rational function of D, was reconstructed from this ex- 
pansion. 

The major part of the calculations was performed on 
the Silicon Graphics Altix 3700 computer (32 Itanium-2 
1.3 GHz processors) using the parallel version of FORM 
US El 13 and took about 18 months in total. The 
diagrams were generated with the help of QGRAF EH- 

It is convenient to introduce the Adler function as 



D(Q 2 



Q 2 d n(Q 2 
6 dQ 2 Q 2 



Q 2 R(s)ds 
(s + Q 2 ) 2 



(4) 



where we have set the normalization scale = Q 2 or 
to /i 2 = s for the Euclidian and Minkowskian functions 
respectively. The specific form chosen in eq. (@J) is par- 
ticularly convenient for perturbative consideration. The 
results for generic values of y, can be easily recovered 
with the standard RG techniques. Since the functions D 
is known to order a^. from Q we only display the O(a^) 
results. 
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The 7ij and n 2 ? terms are in agreement with |2Gj and 

respectively. Once the constants d\ to di are known it is 
a^ matter of straightforward analytic continuation to find 
R, given below for brevity in the numerical form only: 

R = l + 5.6667a s +[35.94- 1.359n/]a 2 

+ a 3 [164.14 - 25.77 n f + 0.259 n)} (6) 
+ aj [39.34 - 220.9 n f + 9.685 n) - 0.0205 nj] . 

In order to better understand the structure of the a* 



term in © it is instructive to separate the genuine five- 
loop contributions from the Adler function D from essen- 
tially "kinematical" , so-called it 2 -terms originating from 
the analytic continuation. We have in mind that for a 
given order in a s these extra contributions are completely 
predictable from the standard evolution equations ap- 
plied to the "more leading" terms in D proportional to 
some smaller powers of a s . The corresponding expression 
for R reads 



R = l+5.667a g +a 2 [51.57 - 15.63 - Tt/ (1.907 - 0.548 )] +af [648.7 - 484.6 - n/(63.74 - 37J7) + n) (0.929 - 0.67)] 
+ a! [9470.8 - 9431.4 - n f (1454.3 - 1233.4 ) + ^(54.78 - 45.10 ) - n^O.454 - 0.433 )] , (7) 

I 



where we have underlined the contributions coming from 
analytic continuation. 

The inclusion of the 7r 2 -terms from higher orders thus 
does not necessarily improve the quality of the perturba- 
tive prediction for the scalar correlator. It remains to be 



seen whether a similar pattern arises in the case of the 
contour improved perturbation theory |2ll |22| applied to 
the r-lepton decay rate. 

At last, specifying nt = 5 we get the corresponding 
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generalization of eq. 10 

R = 1 + 5.6668 a s + 29.147 a 2 + 41.758 a 3 s - 825.7 a A s 
= 1 + 0.2075 + 0.0391 + 0.0020- 0.00148. (8) 

In the last equation we have substituted a s — 0.0366 
which corresponds the Higgs mass value Mjj = 120 GeV. 
The comparable sizes of the C( a s) and the 0(a 4 ) terms 
(the fourth and the fifth terms in eq. ©) may be inter- 
preted as a consequence of the accidentally small coeffi- 
cient for the term. 

COMPARISONS TO PREVIOUS ESTIMATIONS 

Following [23, predictions for higher order coefficients 
of QCD correlators are usually made first for an Eu- 
clidian quantity and the result for the corresponding 
Minkowskian one is then derived as a direct consequence. 
On the other hand, a quite remarkable feature of our re- 
sult for R are almost complete cancellations which take 
place for every power of nj between the genuine five- loop 
contributions and the "trivial" ones from analytic contin- 
uation as illustrated by the corresponding decomposition 
displayed in eq. Q . This fact alone might indicate prob- 
lems for the traditional way of obtaining predictions for 
R. With the exact result in hand one can easily check 
this suspicion. 

Indeed, in Table Q our results are compared with pre- 
dictions obtained in works [3 HH IHl ■ Note that the 
Principle of Minimal Sensitivity (PMS) |27j| and that of 
Fastest Apparent Convergence (FAC) |28j| used in [24| 
produce identical result at order a 4 . Also note that 
the Asymptotic Pade-Approximant Method (APAM) was 
utilized in 25] to produce the prediction directly (and 
only) for the absorptive part R. In contrast, the method 
of Naive NonAbelianization (NNA) [2^| has been applied 
in 26] to the Adler function only. The two predictions of 
FAC/PMS for P4 correspond to either the consequence 
of the prediction for dn (the fifth line) or to the direct 
application of FAC/PMS to estimate (the sixth line). 
As a consequence of the large cancellations in f 4 the sec- 
ond prediction looks much better than the first, despite 



the fact that the estimation of the Euclidian coefficient 
e?4 is quite close (within 10%) to the exact result. In 
fact, at order aj? the cancellations in questions are much 
less pronounced with the result that the corresponding 
prediction for r^, obtained from d^, is significantly more 
accurate than the direct estimation of f 3 . NNA pre- 
dictions have correct signs and sensible magnitudes as 
observed earlier j^. Finally, the APAM estimation of 
T4 fails to reproduce even the sign of the exact result. 

Predictons of the prescripti on p roposed by Brodsky, 
Mackenzie and Lepage (BLM) |30| for the n/ dependent 
terms of order a 4 have been communicated to the au- 
thors 36]: a*(-260n/ + I3n 2 f - 0.046nJ) and are also 
in reasonable agreement with the exact result of eq. 



TABLE I: Comparison of the results obtained in this paper 
with earlier estimates based on PMS, FAC and APAM. 



n f 


3 


4 


5 


di (exact) 


5588.7 


4501.1 


3512.2 


d± ([24], PMS. LAC) 


5180.3 


4093.0 


3100.5 


<L ([261. NNA) 


1592.8 


1521.4 


1484.1 


fi (exact) 


-536.8 


-690.7 


-825.7 


h ([24], PMS, LAC) 


-945.2 


-1098.8 


-1237.4 


h ([24], PMS, FAC) 


-527.8 


-748.6 


-949.4 


fi (1251. APAM) 


252 


147 


64.2 



QUARK MASS BOUNDS 

As an application of our result for the scalar correlator 
we consider the well-known bounds for the light quark 
masses 0, . The constraints follow from the known 
values of the it or K pole contributions and the positivity 
of the spectral function R and depend on the scale Q 
used in evaluation of the polarization operator. We are 
going to discuss two types of bounds suggested in [32| . 
the "quadratic" bound 



[m s (Q) + m u (Q)Y 



> 



16tt 2 2/fM 4 




2^ (Q 2 )-|(l + M|) 


!^i(0 2 ) + ^l 


(l + #) 


iV 2 (Q 2 ) 


N c Q 4 , 




, 5 3^ (Q 2 )^2(Q 2 )-2(^(Q 2 )) 2 



(9) 



and the "linear" bound Here the functions and T2 are defined as (normal- 

_, ized to one at the leading order, that is cq = 1/6, ci = 

K(Q) + m , ( ^>lg ^l 1+ j) . (10) 
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-1/6, c 2 - 1/12) 



Fn(Q ) = Cn (Q ) 



2\l+n 



dQ 5 



2+n 



n(Q 2 ). (ii) 



The results of both eqs. I|9I1U|I for m s + m u depend on 
the choice of Q and can be transformed to bounds for 
(m s +TO„)(/i) for any (not too small) \x using the standard 
mass evolution equation. Following [32l | we will use fj, = 
2 GeV as the reference point. 

Let us start from ©. Using as reference value 
a s {M T ) = 0.334 and the standard evolution equations 
for the mass and coupling constant we get (the indices 4 
and 3 stand for the order in a s ) 

[(m s + m u )% ( ^ GeV^j > (103 MeV) 2 , (12) 

which should be compared to the three-loop bound [13] 



[(» 



m u ) 2 } 3 



H=2 
Q=2 



GeV > (111 MeV)' 



(13) 



The bound l|13l) . if valid, is already in conflict with 
significantly lower value 75(8) MeV derived by one of 
lattice collaborations (3j|. (However, for a significa ntly 
larger lattice result around a hundred of MeV, see |3J, 

E3). 

The problem was "solved" in [6j by observing that the 
quadratic combination of the T% functions appearing in 
the denominator of © displays a poor pattern of conver- 
gence: 



3^0-^2 - 2(^i) 2 = 1 + 0.83 + 0.61 + 0.51 



(14) 



and then suggesting to increase Q in (|13|l till, say, 
2.5 GeV, which lowers the rhs of JT3|) to 86 MeV and 
the rhs of to 81 MeV. 

From our point of view it is more relevant to consider 
the convergence pattern of the whole ratio appearing in 
the rhs of eq. Q- Indeed, after expanding 10 in a s and 
disregarding all terms in a s higher than we get 



[{m s + m u )\ ( GeVj > (179 MeV) 2 x 
{ 1 - 6.44 a s - 12.83 a 2 + (482.95 - 525.2 ) a\ 



(-2561.8 + 6948.34 - 



4439.9 



(15) 



In order to demonstrate the importance of higher order 
corrections we have underlined in H15(l all terms originat- 
ing from the contributions of order a 3 (eventually multi- 
plied by a subleading 0(a s ) term) to the functions J-q,T\ 
and T% and, correspondingly, boxed the one originating 
exclusively from O(af) terms in the same functions. Af- 
ter substituting the value for a s (2 GeV) = 0.312 into 
eq. I|15f) we arrive at 



to be compared to 

[{m s + m„)] 3 (/i = 2 GeV) > 79 MeV. 



(17) 



If, on the other hand, we would try to use eq. (|15fl with 
only the boxed term nullified we would get an astonish- 
ingly different value 

[(m s + m u )] 3 (fi = 2 GeV) > 141 MeV. (18) 

In similar way one can examine the second bound i|10[l . 
For the lowest choice Q = 1.4 GeV used in [3^ we im- 
mediately get 

[(m s + m u )]{n = 2 GeV) > 76 MeV and 78 MeV, 

with the numbers corresponding to the use of five- and 
four-loop expressions for the function Tq correspondingly. 

On the other hand, the expanded version of eq. i|10|) 
reads: 



■ m„) 2 ] 4 



fi=2 
0=1.4 

. 2 



GeVj > (106 MeV) 2 x 
{ 1 - 3.67 a s - 0.73 a 2 s + (54.7 - 77T) a 3 s 



511.8 



(19) 



+ (-190.1 + 567.4 
It leads to the following bounds 

[(m s + m u )] 4 (p = 2 GeV) > 72 MeV, (20) 



[(m s + m„)] 3 (/i = 2 GeV) > 74 MeV. 



(21) 



[(' 



m u )] 4 (n = 2 GeV) > 77 MeV 



(16) 



Again, if we were using use eq. I|19f) with only the boxed 
term nullified we would arrive at significantly different 
value 

[(m s + m u )] 4 (/x = 2 GeV) > 81 MeV. 



CONCLUSION 

We have computed the correction of order to the 
correlator of quark scalar currents in the massless limit. 
Our result demonstrates a remarkable interplay between 
the genuine five-loop terms and the effects due to the an- 
alytical continuation. The newly computed contribution 
stabilizes the (quadratic) quark mass bound of |32j and 
pushes it significantly down, thus, avoiding any potential 
conflict with lattice results. 
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